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Disturbances in a Multi-Velocity Plasma’ 


J. R. PIERCE ann J. A. MORRISONT 


Summary—The single-velocity treatment of ac disturbances in 
electron flow is often accurate enough even when the flow actually 
has an infinitely broad but peaked velocity distribution, such as 
a Maxwellian or a bell-shaped distribution. In seeking a better 
approximation, one is tempted to expand in terms of the moments 
of the velocity distribution and to disregard terms beyond that 
involving the second moment. This leads to a dispersion equation 
describing waves which strictly do not exist. A linearized 
analysis by transform methods predicts a field which oscillates in 
an exponentially damped manner. This is not the oscillation of a 
normal mode, however, for various velocity classes of the charge 
distribution have ac densities which grow with time. This indicates 
that the linear expression will hold over a finite interval only. 

A simple example which compares a multi-stream analysis with 
an analysis based on the charge density in phase space (the distri- 
bution function) indicates the same sort of failure of both of the 
linearized theories at large times. Thus, it appears that the failure 
is characteristic of the problem rather than of the method of solu- 
tion. The method of solution is a matter of choice. 


INTRODUCTION 


ORKERS in the field of electron tubes have had 
a number of years to deal with problems con- 
cerning the multivelocity nature of actual 
electron flow. Some of these problems are concerned with 
noise in electron beams. Here the crucial effect of multi- 
velocity flow lies in the fact that the spread in electron 
velocities, together with the discrete nature of electrons, 
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leads to two independent sources of noise near the cathode, 
and these noise sources give rise to noise space-charge 
waves in the electron beam.*” 

In noise calculations, one can sometimes get sufficiently 
accurate results by disregarding the effect of the multi- 
velocity nature of the electron flow in dealing with the 
propagation of noise disturbances along the electron 
beam, but this can lead to serious errors in the low- 
velocity region near the cathode? Thus, we are faced 
with the questions: How can we best take into account 
the actual velocity distribution? Is there any simple 
approximation which will help us? It is the purpose of 
this paper to review our present knowledge of this matter. 
In the course of this review, we shall exhibit the general 
features of the solution of multi-velocity problems by 
both the multi-stream approach and the distribution- 
function or density-function approach.* 

We shall consider the case of disturbances in an infinite 
cloud of electrons; we may imagine their average charge 
to be neutralized by a fixed positive charge. We shall 
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consider disturbances in which the fields vary in the z 
direction only. The treatment will be Newtonian and 
electrostatic. 


Waves In Mvuuri-Vetoctry Flow? 


Let us consider an undisturbed cloud of electrons in 
which the average charge density associated with electrons 
whose z-component of velocity lies in the range du at % is 
—dpo. Let —p, be the total charge density for electrons 
of all velocities. Then we can write 


dpo = pou) du. D 


Here P(u) describes the velocity distribution of the 
electron charge density. It is clear that the form of the 
expression has been chosen so that 


[ Pù du = 1. 9) 

How shall we take a disturbance in such an electron 
cloud into account? Let us consider a case in which all ac 
quantities vary with time and position as 


eiet, 


It is easily shown that for electrons which have a velocity 
u in the absence of the disturbance, the ac convection current 
di associated with the field Æ is given by 


= jole/mE dpo, 


ae (w — Bu)’ 


(3) 
To get the total convection current 7 we must integrate 
over all velocities; thus 

; ? Plu) 

jwle/ m) pE d lw a Bu)’ du. 

A plausible approximation has been proposed.’ Let us 
assume that @u is very small compared with w, and expsnd 
the denominator in the integral so that 


i= 


(4) 


es ZIE m) pol? 


[+ 2Gu/s) + B(Gu/e)* + Pi) du 
j= SEE 1 4 Aau lu) + 306/0)%U") + --). (8) 


This procedure is strictly justified if |@u/w! is less than 
unity for the greatest value of |u| in the velocity dis- 
tribution. However, it is dubious to use this expansion 
in connection with infinite velocity distributions. 

If we assume (u) = Q (zero average velocity) and 
disregard all terms beyond the term in (u*), as we 
might if |6u/w| were very small for the largest |u| in the 
distribution, we obtain 

i = HEME + 306/02). 


w 


(6) 
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Now let us assume that the total electric field is due 
to the convection current 7. As the total current density 
must be zero 


t + josh = 0. (7) 
From (6) and (7) we obtain 
w = o + 3(8/«)*(u")). (8) 
The quantity wp, called the plasma frequency, is 
op = EAD gy, 9 
Eq. (8) can be rewritten 
( he T a 2 1/2 
B= = 302) ((2) — 1) ; (10) 


This relation is usually known as the dispersion relation 
of Bohm and Gross.” 

It appears that a disturbance of a particular phase 
constant (and wavelength) is associated with an exci- 
tation of any frequency w greater than the plasma fre- 
quency. If the mean square velocity (u) approaches 


_ zero, however, the range of frequencies corresponding to 


any finite range of phase constants 8 approaches zero at 
= Wy. 

To explore this matter, let us assume a particular 
velocity distribution which is simple enough to handle 
analytically. Let 


Ludu, 
mu” + uo) 
Here —p is the average charge density for all velocity 
classes of electrons. Thus, from (3) and (11), the total 


convection current will be 


i = —jle/m)(poiu/T)0/ BE 


(11) 


dpo 


; i du = 
dae (UU + jito) ~ juotu — w/B)? 


The integral can be evaluated by contour integration. 
As the integrand goes to zero at infinity in all directions 
in the complex u plane, we can close the contour either 
above or below the real axis. If we assume w to have a 
positive imaginary part, it is most convenient to close the 
contour below the real axis, including only the pole at 


(12) 


u = —juy. (13) 
This leads to 
a —jle/m) powt 
o (w + jB a4) 


If, however, we assume w to have a negative imaginary 
part, we obtain 
_ Ziem peek 
em 78%)” 
At this point, we look for a solution of both (7) and 
(14) or (15). From (7) and (14) we obtain 


(15) 


(16) 


The imaginary part of w as given by (16) is negative, 
contradicting the assumption which led to (14), that the 
imagmary part of w was positive. We arrive at a like 
contradiction if we use (15) and (7). Thus, there are no 
normal modes in which all the ac quantities vary as exp 
(Got — 782). 

We might think that we fail to find waves of the sort 
described by (10) because for the velocity distribution 
(11), <u’) is infinite. If we believed relation (8), this 
would. lead to an infinite frequency of oscillation for an 
oscillation of finite wavelength (8 = 0). However, we shall 
show that despite the infinite mean square velocity, the 
distribution exhibits a calculable and reasonable behavior 
following excitation. Further, this behavior is similar to 
that which Landau’ has found for a plasma having 
a Maxwellian velocity distribution. Also, Walker’ has 
studied the general case of waves in streams with an 
arbitrary velocity distribution. He finds that when the 
velocity distribution has a finite width, there are un- 
attenuated waves having phase velocities lying outside 
the range of electron velocities. When the velocity dis- 
tribution has an infinite width, there are no inattenuated 
waves; there may in this case be growing and decaying 
waves if the velocity distribution has two or more peaks 
with sufficiently deep valleys between them. 

Thus the dispersion relation (10), although plausible, 
has no strict significance in the case of infinite velocity 
distributions, although it does give the correct phase 
constant for small values of (8/w) (w) in the case of 
the Maxwellian velocity distribution.” We may ask, 
however, what sort of disturbances can exist in a plasma 
with an infinite, single-peaked velocity distribution? 

We are unable to treat this problem in its true nonlinear 
form. It is possible to obtain a linearized solution, pre- 
sumably valid for small disturbances. We will find, 
however, that linearized solutions of electron-flow problems 
can be in error at large enough times even when the initial 
amplitude of excitation is small. 


o= zhw, =< Buo. 


A PLASMA OF [INTERACTING PARTICLES 


Let us consider what happens when a plasma of particles 
with the velocity distribution of (11) is acted on by an 
impulsive force with a sinusoidal spatial distribution and 
by space-charge forces. For convenience, we will write the 
initial force on a particle as —eF, where we may think of 
F as a nonelectrical quantity which is introduced merely 
to impart an initial velocity to the electrons. We will 
assume that F has the form 


—V, 


= (¢/m) 


OF (17) 
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where ô(t) is a delta function. This force will give a particle 
at = 0 an additional velocity v, at £ = 0. Considering the 
stream corresponding to an initial velocity u, we have, 
according to the linearized theory, 


aù 


ð — 
at (dp) = az? 


(di) = (dou — po, 
(18) 


where (dp), (di) and v are the perturbation density, 
convection current, and particle velocity, respectively. 
The last. equation of the group asserts that the total 
current density is zero. If we let do = dp*e™***, etc., (18) 
becomes, upon taking Laplace transforms, 


as A ers -£ Ee 
(S — jeu o = vn — 7 Bs, 


TE + SeE* = 0, 
S(dp*) = joldi”), 


(di) = (dp*)u = (dpo)v*, 


(19) 


where S is the transform parameter. Eliminating v*, 
(dp*) and H*, we obtain 


_ [ s, + oF 
(di*) = —(dpp) “TS jew? 


Inverting (20), and using the convolution theorem,® 


(20) 


di* = 


— (dpo) (a + judi" 


£ m = iBu(t—7) 
-+ m t(D re ar), (21) 


Integrating (21) with respect to u from — œ to + œ, 
where dpo is given by (11), 


it = 


=at — pude 


En a i TAME — De Bue ar), (22) 


This is an integral equation for 7* which may be solved 
by taking Laplace transforms, 


ge a 
-a| Se pi oF] i i S 
(SFe 7 


mE (SE Buy Far? 2) 


8 W, Magnus and F. Oberhettinger, “Formulas and Theorems 
for the Functions of Mathematical Physics,” Chelsea Publishing 
Co., New York, N. Y., ch. 8, sec. 2; 1954. 
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where w, is the plasma frequency as defined in (9). From 
the second equation; of (19), it then follows that 


TF Pop. 1 = 
e 1S + Buy Fall ee 
Inverting (24) we find 
E* = ot CE (25) 


pE 


This bears not the slightest resemblance to the result of 
Bohm and Gross’ (8). Instead of an oscillation of frequency 
other than w, we have a damped oscillation of frequency 
Wp. 

Clearly, although the field has only this exponentially 
decaying sinusoidally oscillating component, currents 
carried by streams of various velocities must have a riore 
complicated behavior. 

Let us then consider the charge density dp associated 
with the electrons of charge density dp, and of ave‘age 
velocity u in the range du. There will be two compon ənts 
of the charge density dp associated with these electrons, 
that due to the impulsive force F, and that due to the 
subsequent electric field Æ, which is zero fort < 0. From 


(19) and (20), 
-F 
m 


— ju 


The first component of (dp*) is that caused by the im- 
pulsive force F att = 0, 


(d) = (dpi) + (dos): 20) 


ev (dpo) s 
d o ‘ 2 
Ge) = g 27) 
Inverting (27), we obtain 
dp, = —jBv,tdp)e™'e"**. (28) 
From equations (26) and (24), 
n jB(dpo) £ E* 
(dp*) = me 
= S — jpu’ 
JBr,(dpow, (29) 


= (S = jw CS + Bu) + oF] 


Expressing the right-hand side of (29) in partial fracticns 
and inverting, we find the part of the charge density due 
to the field Æ is 


2 = ji ? d 0 ae oz 7 -| i 
Na Boe oe — [bu — j 


Twpt 


ca aa 
3 Vo Ap Bu,€ lw, — Bu — jul 


ei eet 


= erl” 
[wp + Blu — juto)] 


(u — juge” 


+ 2 d o 2 26 ibut - - — 
OPPO Cat, = ja) = wh 


(3)) 
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From (28) and (80) we see that there is a component 
of charge density dp, associated with electrons in the 
velocity range du at u which increases linearly with time; 
it is 


_ iB: iBut _ [Bu i juyi? ] 5 

dpa iaj j dpo Gv,le F == [Blu S jo) | . 81) 
We note that (31) says that when 
$ y 2 

_ Begt Bu + Juo) | = ii (32) 


— [Bu — jud | 
the peak value of the ac charge density dp, will be equal 
to the de charge density dp, of the electrons of the stream 
with which it is associated, Clearly, the approximation of 
linearity will have failed before this time. 

Let us consider the case in which the wavelength of the 
disturbance is very short, so that 


| Blu + iud | > wy. (33) 
If this is so, (25) will hold only when 
Bvt K1. (34) 


Thus, for this case, (25) will hold only until the oscillation 
has decayed to a point such that 


ge Fue! sy g tele (35) 


The smaller the initial velocity of excitation », compared 
with the breadth of the velocity distribution, the further 
the oscillation decays in amplitude before its behavior 
is seriously affected by nonlinearities. 

The number of cycles of oscillation which take place 
before nonlinearities become important is also of interest. 
When £ is large and (38) holds, this number Ñ is 


ee ta (36) 


T = P 
A Qn ~~ rhv, 


As a practical matter, nonlinearities may not be par- 
ticularly important. For a plasma with a velocity spread 
comparable to the Maxwellian spread for common 
temperatures, the rate of exponential decay will be small, 
and for times of interest the behavior will be close to that 
of a single-velocity plasma. Further, there is every reason 
to believe that in the exact solution, as in the linearized 
solution, the field will become very small at a sufficiently 
great time subsequent to excitation. 

It is, however, interesting to note that the exponentially 
decaying solution of the linearized problem is not a 
normal mode of oscillation in which all quantities vary 
exponentially with time. In this sense, it is not a wave. 
Initially, the field varies exponentially with time with a 
complex exponent, but the charge densities of various 
velocity classes of electrons do not. Further, it appears 
that ultimately even the exponential variation of the 
field must break down. 

It might be objected that the failure of linearity im- 
plied by (82) is a product. of the nature of the analysis 
rather than of the actual physical nature of the problem. 
The formulation used above is a multi-stream formulation, 
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in which we group together electrons which have the same 
velocity prior to the introduction of the disturbance. It is 
well known that, in such an analysis, difficulties arise 
when the ac velocities cause electrons to overtake one 
another. 

We can avoid this difficulty by working with the charge 
density in phase space (the distribution function) rather 
than with the charge density in configuration space. 

Liouville’s theorem’ states that the total time derivative 
of the distribution function in the vicinity of a given 
particle is zero. By expressing the total derivative as a 
sum of partials, and by expressing the time derivative of 
the absolute particle velocity w by means of the equation 
for the particle acceleration, we find that the distribution 
function f(t, z, w) satisfies the equation 

Lp wle fw emo. (37) 

In dealing with this equation, w is to be regarded as 
one of the coordinates, just like ¢ and z. There is no varia- 
tion of w with ¢ or z; instead, f varies with ż, z, and w. 

It is assumed that at t = 0 we have f = P(w), a func- 
tion of w alone, where P satisfies (1) and (2). We linearize 
(87) by replacing df/dw by dP/dw. Let g(t, w) be defined by 


f = Pw) + e g(t, w). (38) 


With the substitution of (38), where F is as given by (17), 
the linearized equation becomes 


T E hee, JE 
20 ipug = | Ln» aaO |S (39) 
where E = E*e 1%, Taking Laplace transforms, 
: dP 
eds Fi = * — . 
(S — j8w)g E E — Jg To (40) 
Also we have the equations, 
ôe ot, OE 
Ja a ey ae ay 
Letting p = p*e i", i = 01" and taking Laplace trans- 
forms, 
Sp* = jbi*;  i* + SeE* = 0. (42) 
Eliminating E* and 7* from (40) and (42), 
e 
pa hm pe dP, 
(S — jaw) dw a 
Inverting this equation, 
= éj * jRwlt—7r) =x a dP 
g (£ Be Í p*(re dr — v,e Ti (44) 


But, from (88) and the definition of the distribution 
function, 


? A derivation of Liouville’s theorem direct from Newton’s laws, 
without the use of Hamilton’s equations, is to be found in J. 
A ete and Design of Electron Beams,” D. Van Nostrand 

0., Sec. 


Pierce and Morrison: Disturbances in a Multi-Velocity Plasma 


235 


p* = -o f g dw. 

Substituting for g from (44) into (45) and performing 
the integration with respect to w, we have for the bell- 
shaped distribution, 


(45) 


t 
j= (e; 1 pe(a)(t — e dr + jnst"), 
0 
(46) 
This is an integral equation for p* which we may solve by 
taking Laplace transforms, 


eR = (w2p* + JB p00) Tk el his JB rPo (47) 


eee St Bu) P [CSF buo E o] 


From (42) and (47), we derive the same expression for 
E* as was found in (24) by the multistream approach. Let 
us now determine the distribution function perturbation. 
From (48) and (47), 


oe __ & — OP 
I= “(S = jw) dw 
2, dP 
Wolo Tiy 


gtg (48) 


+ ONS + ha Fel = 
The first component g, is that caused by the impulsive 
force F at t = 0. We obtain it from (40) by disregarding 
the first term in the brackets on the right and inverting 
for g. This gives 


g = — 0,0"! (49) 
From (48) 
2 aP 
_ Wty Te 
o = jaws + Buy F ol] a 
„dP 
. Oo Tu ( 1 (S++ j8w + 26s) 
[B luo + jw)? + ws] S — jBw) KS + Bur)” + wl 
Inverting, 
dP 
Ynn dw iwit —fuot 
go = kope- mi lwp t — & wp COS (wpt) 


-+ Bluo + jw) sin @,t))). 


Now in linearizing the distribution equation, we replaced 
ðf/ðw by dP/dw in (89). Since f is given by (88), the 
linearization is valid only if 


(51) 


<8 X | æl | t>0 andal w. (52) 
From (49) and (51) it follows that 
m a — ie | dP | 
ra ea <| aw | 3) 
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This implies that both 


vB | w — juo |? EP | | dP | o 
| Ws 7 B(w ag jo)” | dw? | < dw l? all Ww, 54) 
a restriction on the initial distribution function, and 
x 2 


| op — B(w — juo) | 

The latter condition is identical with that obtained for 
the validity of the linearization in the multi-str3am 
approach. Thus, the limitation imposed by nonlinearity 
is seen to be inherent in the problem, and it does not 
represent a weakness of the multi-stream methoc of 
solution. 

We finally discuss an alternate approach to the trans- 
form solution. In the distribution function approach, 
(43) was first inverted to give g; this value of g was then 
substituted into (45) and the integration with respect to 
w carried out. This led to an integral equation for the 
spatial density which was then solved by a further appli- 
cation of Laplace transforms. Alternatively, we could (ake 
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the transform of (45) to find 
p* = — Po I g dw 


[ei _, |f 1 a 
T a| Be nl gt, 


from (43). For the bell-shaped distribution, it is found by 
contour integration that 


(56) 


r 


-j8 
P i ga (ES OA 
l. (S — jw) dw z _ 5g 
(S = Bu)?’ R(S) < 0, 


where R denotes the real part. Thus it appears that the 
change of order of integration, with respect to w and with 
respect to ¢, in obtaining p* from (45) is not valid for all S, 
From the nature of the Laplace transform and the in- 
version theorem,* the correct transform of p* is given with 
R(S) > 0 and from (56) and (57) we then arrive at (47). 
A similar procedure may be adopted in the multi-stream 
approach. Thus (20) may be first integrated with respect 
to u from — œ to + , and (23) is obtained for RUS) > 0. 


